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Abstract 

We investigate generalizations of the Charlier polynomials on the lat- 
tice N, on the shifted lattice N+1-/3 and on the bi-lattice NU (N+1-/3). 
We show that the coefficients of the three-term recurrence relation for 
the orthogonal polynomials are related to solutions of the fifth Painleve 
equation Pv (which can be transformed to the third Painleve equation). 
Initial conditions for different lattices can be transformed to the classical 
solutions of Pv with special values of the parameters. 
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1 Introduction 

In this paper we investigate the recurrence coefficients of orthogonal polyno- 
mials, in particular generalized Charlier polynomials, and show that they are 
related to the solutions of a continuous Painleve equation when viewed as func- 
tions of one of the parameters. We are interested to see how the properties of 
the orthogonal polynomials are related to properties of transformations of the 
Painleve equation. 

The paper is organized as follows. In the introduction we shall first review 
orthogonal polynomials for a generalized Charlier weight on different lattices 
and their main properties, following [10]. Next we shall briefly recall the fifth 
Painleve equation and its Backlund transformation. We then derive the Toda 
system and show that the recurrence coefficients can be expressed in terms of the 
solutions of the fifth Painleve equation. Finally we study the initial conditions 
of the recurrence coefficients for different lattices. 

1.1 Orthogonal polynomials for the generalized Charlier 
weight 

One of the most important properties of orthogonal polynomials is a three-term 
recurrence relation. For a sequence (p n (x)) n eN of orthonormal polynomials for 



a weight w on the lattice N = {0, 1, 2, 3, . . .} 

oo 

^2p n {k)pk(k)w(k) = 6„ tk , (1) 

fe=0 

where S n ^ is the Kroneckcr delta, this relation takes the following form: 

xp n (x) = a n+1 p n+1 (x) + (x). (2) 

For the monic polynomials P n related to orthonormal polynomials p n (x) = 
7„x n + • • ■ with 

oo 

= Y. P n(k) W (k), 

7 ™ fc=0 

the recurrence relation is given by 

xP n {x) = P n+ i(x) + b n P n (x) + a^P„_i(x). 

The (classical) Charlier polynomials ([4, Chapter VI], [10]) are orthogonal 
on the lattice N with respect to the Poisson distribution: 

°° a k 

^ C n (k; a)C m (k; a)— = a~ n cxp(a) n\ 5 n>m , a > 0. 

fe=0 

The weight Wk = w(k) = a k /kl satisfies the Pearson equation 

Vw(x) = ^1 J w(x), 

where V is the backward difference operator 

V/(ar) = f{x) - f{x - 1). 

Here the function w(x) = a x /T(x + 1) gives the weights Wk = w(k) = a k /k\. 
The Pearson equation for the Charlier polynomials is, hence, of the form 

V[er(x)u;(x)] = t(x)w{x) (3) 

with cr=l and r a polynomial of degree 1 . 

Recall that the classical orthogonal polynomials are characterized by the 
Pearson equation (3) with a a polynomial of degree at most 2 and r a poly- 
nomial of degree 1. Note that in (3) the operator V is used for orthogonal 
polynomials on the lattice and it is replaced by differentiation in case of orthog- 
onal polynomials on an interval of the real line. The Pearson equation plays an 
important role in the theory of classical orthogonal polynomials since it allows 
to find many useful properties of these polynomials. It is known that the recur- 
rence coefficients of the Charlier polynomials are given explicitly by a„ = na 
and b n = n + a for n £ N. 
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The Charlier weight can be generalized [10]. One can use the weight function 

_ r(/3K 

W[X> T(p + x)T(x + l) 

which gives the weight 

k 

w k = w{k) = a > 0, (4) 

on the lattice N. The Pearson equation is now of the form 

_ , , a — x(B — 1) — x 2 , . 

Vw(x) = — '- w{x). 

a 

The monic orthogonal polynomials P n (x;a,f3) for the weight (4) satisfy 

00 k 

Y,Pn(k;a,(3)P m (k;a,(3)-^- =0, n ^ m, (5) 
k=o [P)kK - 

and the three-term recurrence relation is 

xP n (x;a,(3) = P n+1 (x; a, f3) + b n P n (x; a, (3) + a^P n _i(a;; a, /3), 

with initial conditions P-\ = and Pq = 1. If we replace a by /3a and then take 
the limit (3 — > oo, one recovers the Charlier polynomials [10]. 

Theorem 1.1. [10, Th. 2.1] The recurrence coefficients for the orthogonal 
polynomials defined by (5) for the weight (4) on the lattice N satisfy 



b n + b n - 1 -n + /3= (6) 
at 



,2 „\/„2 



(a n+1 - a)(a n - a) = a(b n - n)(b n - n + [3 - 1) (7) 
with initial conditions 

, V"^ 2 ^) /ON 

where Ip is the modified Bessel function. 

The system (6)-(7) can be identified as a limiting case of an asymmetric 
discrete Painlevc equation. In this paper we show that it can be obtained from 
the Backhmd transformation of the fifth Painleve equation. 

One can use the weight (4) on the shifted lattice N + 1 — j3 and one can also 
combine both lattices to obtain the bi-lattice NU (N+ 1 — f3). The orthogonality 
measure for the bi-lattice is a linear combination of the measures on N and 
N+l-fi. The weight w in (4) on the shifted lattice N+l-/3 = {l-/3,2- 
(3, 3 -/?,...} is given by 
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We see that up to a constant factor it is equal to the weight on the original lattice 
N with different parameters. The corresponding monic orthogonal polynomials 
Q n on the shifted lattice satisfy 

oo 

Qn(k + 1 - P)Q m (k + 1 - P)v k = 0, n ^ m, 

fc=0 

and one has 

Qn(x) = P„(x + /3-l;a,2-/3). 

In [10, Th. 2.2] it is shown that the recurrence coefficients in the three-term 
recurrence relation 

xQ n {x) = Q n+ i{x) + b n Q n (x) + a^Q n -i{x) 

satisfy the same system (6)-(7) (with hats) but with initial conditions 

a - U, b - — . (10) 

ll-/3(2y/a) 

Using the orthogonality measure /i = c\\i\ + C2/A2, where ci, ci > 0, [i\ 
is the discrete measure on N with weights = w(k) and H2 is the discrete 
measure on N + 1 — /3 with weights = w(k + 1 — /3), one can study the monic 
orthogonal polynomials R n (x) — R n (x; a, ft, r), where r = Cq,/c\ > 0, satisfying 
the three-term recurrence relation 

xR„(x) = R n +i(x) + b n R n (x) + a 2 n R n -i{x). 

The orthogonality relation is given by 

00 00 

Ci ^ Rn(k)R m (k)w k + C 2 53 + 1 - P)Rm(k + 1 - P)l> k = 0, 771 ^ 71. 

According to [10, Th. 2.3] the recurrence coefficients and b n satisfy the system 
(6)-(7) (with tilde) but with initial conditions 

a = 0, b = y/a- fe) r , , T 77T1=\- ( n ) 

lp-i[2y/a) + Tl 1 -p(2y/a) 

Thus, it is shown in [10] that the orthogonal polynomials for the generalized 
Charlicr weight on the lattice N, on the shifted lattice N + 1 — (3 and on the 
bi-lattice N U (N + 1 — 0) have recurrence coefficients a\ and b n which satisfy 
the same nonlinear system of discrete (recurrence) equations but the initial 
conditions are different for each case. 
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1.2 The fifth Painleve equation and its Backlund trans- 
formations 

The Painleve equations possess the so-called Painleve property: the only mov- 
able singularities of the solutions are poles [7]. They are often referred to as 
nonlinear special functions and have numerous applications in mathematics and 
mathematical physics. 

The fifth Painleve equation Py is given by 

1 1 y' Av- 1 ) 2 ( a.. , B \ , °y , D v(y + 1 ) 



where y = y(t) and A, B, C, D are arbitrary complex parameters. 

The fifth Painleve equation Py in case C^O, D = can be reduced to the 
third Painleve equation [7]. There exists a Backlund transformation between 
solutions of the fifth Painleve equation with D = 0. 

Theorem 1.2. [2], [ff] Let y = y(t) be a solution of Py with parameters 
A, B, C, D given by B = -p 2 /2 / 0, C ^ and D = 0. Then the function 
2/i = Vi (*) given by 

2Cty 2 (y-l) 

2/1 = 1- 



2A|f»(y-l)2-(p(y-l)-V) s 



is a solution of Py with parameters A\ = A, B\ = — (p — l) 2 /2, C\ = C, D\ 
0. The inverse transformation is given by 



y = i- 



2Cty 2 { Vl - 1) 



2Ayl{y 1 -iy-{{l-p){y 1 -l)-ty' i r 

In case there exists another well-known transformation. 

Theorem 1.3. [7, Th. 39.1] If y — y(t) is the solution the fifth Painleve 
equation (12) with parameters A, B, C, D, then the transformation 

T si,e 2 ,e 3 ■ 2/ — ^ 2/1 

gives another solution y\ = yi(t) with new values of the parameters A\ , B±, C\, D\, 
where 

2dty 

2/1 = 1 



ty' - ay 2 + (a - b + dt)y + b' 



Al = ~^D {C + d{1 ~ a ~ 5))2 ' Bl = ^D { ° ~ d{l ~ a ~ &))2 ' 
Ci = d(b - a), Di = D, 

a = e 1 V2A, b = £ 2 \ /= 2B, d = e 3 V-2D, e 2 = 1, je {1,2,3}. 
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Example. Let y = y n (t) be a solution of Py with parameters A = —B = 
n 2 /8, C = 0, D = —8. Then the function 

= 1 l&tey 

Vn+1 ~ 2ty' + ny 2 + 8tey - n 

is a solution of P v with A = —B = (n + l) 2 /8, C = 0, D = -8, e 2 = 1. 

Remark 1. For our purposes it is sufficient to use the standard Backlund 
transformations of the Painleve transcendents which are given in N1ST Digital 
Library of Mathematical Functions (DLMF project 1 ), see [9, Chapter 32]. This 
also contains various algebraic aspects of the Painleve equations. The inter- 
ested reader can easily re-formulate the transformations within the framework 
of Noumi-Yamada's birational representation of the afhne Weyl groups [8] . 



1.3 Relation to the third Painleve equation 

The third Painleve equation Pm is given by 

u " = _ + hay? + ft + ju 3 ~ -, (13) 

u z z u 

where u — u(z) and a, j3, 7, S are arbitrary constants. By a scaling of the 
independent variable, the parameters 7 and S can be taken equal to 1 and — 1 
respectively. There is a transformation between solutions of Py and Pm. 

Theorem 1.4. [7, Th. 34.1] Let u = u(z,a,(3) (7 = 1, 5 = —I) be a solution 
of Pm such that R = u' — eu 2 — [as — l)u/z + 1^0. Then the function 

2 

y(t) = 1 e 2 = 1, 

is a solution of Pv with parameters 

A= <fi-te+*r tB= _&+a*-*r tC= _ etD=0 . 

32 32 
The inverse transformation is given by the following theorem. 

Theorem 1.5. [7, Th. 34-3] Let y(t) be a solution of Py with parameters 
A, B, C 2 = 1, D = 0. Then the function 



u(z) = ^y, z 2 = 2t, $(t) = ty - V2Ay 2 + (V2A + V-2B)y - V-2B ± 
is a solution of Pm with parameters 



= 2C{V2A - V^2B - 1), (3 = 2(v / 2A + v /r 2B), 7 = 1, 5 = -1. 



1 http: //dlmf . nist .gov /32 . 7 
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2 Main results 



In this section we relate the recurrence coefficients to the solutions of the fifth 
and third Painlevc equations and identify the discrete system (6)-(7). 

2.1 The Toda system 

It is well-known that for a measure of the form e tx dfi(x) the recurrence coeffi- 
cients satisfy 

K)' : = -77 (On) = a l( b n ~ K-\), , 

d a (14) 

with initial conditions a^(0) and 6„(0) given by the recurrence coefficients of the 
orthogonal polynomials for the measure fi. The weight (4) is of this form when 
we put a = e*. Hence, the recurrence coefficients of the generalized Charlier 
polynomials on the lattice N, when viewed as functions of the parameter a with 
respect to weight (4), satisfy 

al 




{K - b„_i), 
\ a n+l a n>- 



a (is) 



We get this system from (14) by applying the transformation a = exp(t) (see 
[.')]). The weight (9) is also of this form when a = e', since 

w(k + l-p)~ [P> 



r(2-p)k\(2-p) k 

Hence, the recurrence coefficients for the shifted lattice also satisfy the Toda 
system (15). Also for the bi-lattice we have that w(x) = e t2 T(^)/r(/3+a;)r(a;+l) 
when a = e , which is of the required form, so that the recurrence coefficients 
again satisfy the Toda system (15). The only difference between these three 
cases is in the initial conditions which correspond to the recurrence coefficients 
of the generalized Charlier polynomials with n = on respectively the lattice 
N, the shifted lattice N + 1 - /3 and the bi-lattice N U (N + 1 - /3). 

2.2 Relation to the Backlund transformation of the fifth 
Painleve equation 

First, let us obtain a nonlinear discrete equation for b n (a). Denoting a\ = x n , 
we have from equation (7) that 

a(a + n(fi - n - 1) + (1 + 2n - /? - b n )b n - x n ) 

X n +1 = ■ (16) 

Equation (6) gives 



f3-n + b n - 1 +b r , 
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Substituting the last equation and (16) into equation (6) with n + 1 we have a 
nonlinear difference equation of the form 

F(b n -i,b n ,b n+1 ,a) = 0, (17) 

for some F which we do not write out explicitly, since it is quite long. We shall 
show later on that equation (17) can in fact be obtained from the Backhand 
transformation of the fifth Painleve equation. 

To derive the differential equation for b„(a) we use (16) and from equation 
(6) we have 

(XTl | TTiXji ^Xfi ^n'^'Ti 

On-l = — ' 

•En 

Substituting this into the first equation of (15) we get 

, (f3 — n + 2b n )x n — an 
x n — 

a 

and from the second equation of (15) we have 

, a(n — a + n 2 — n[3) - a(l + 2n — (3)b n + ab 2 n + 2ax n — x 2 n 

°n - \ ■ ( 1S J 

a(x n - a) 

Differentiating equation (18) and using the expression for x' n , we get an equation 
for b' n , b n and x n . Finally, we can eliminate x n by computing the resultant 
of this equation and equation (18) and we get a nonlinear second order and 
second degree equation 

G(b';,b' n ,b n ,a) = 0. (19) 

Now the difficulty is in identifying this equation. Let us scale the indepen- 
dent variable b n (a) = B n (a/ki) and take a = kit. The equation (19) becomes 
Gi(-B", B' n , B n , t) = 0, where the differentiation is with respect to t. The fol- 
lowing result holds. 

Theorem 2.1. The equation Gi{B^,B' n ,B n ,t) = is reduced to the fifth 
Painleve equation P\ by the following transformation 

R M l + n+(P-3n-2)y+(l + 2n~ [3)y 2 + ty> 

- My-1) ' (20) 

where y = y(t) satisfies Py with parameters given by 

A-V^f, *~<2±2. C = 2h, D-0. (21) 

Remark 2. It is known [7] that Py with D = can be reduced to the third 
Painleve equation. However, to show that equation (17) can be obtained from 
the Backlund transformation, it is more convenient to use Py. 

Remark 3. The parameters (21) are invariant under fj — » 2 — (3; compare with 
the weight on the shifted lattice (9). 
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Let k\ = 1 and a = t. We use Theorem 1.2 to get y n +i and from 
y = y n . In particular, 

, 4%-l)y 2 

2/n-l = 1 



(P - l) 2 (y - 1)V - ((n + - 1) - V) 2 
is a solution of P v with A = (f3 - l) 2 /2, B = -n 2 /2, (7 = 2, D = and 

4t(y - l)y 2 



Vn+l = 1 - 



(p - l) 2 (y - 1) V - ((n + l)(y - 1) + V) 2 



is a solution of P v with A = (p - l) 2 /2, B = -(n + 2) 2 /2, C = 2, D = 
provided that y = y(t) is a solution of Py with (21). Using (20) we can express 
b n ±i in terms of y' ', y. For instance, 

_ 1 + n + (/3 - 3n - A)y + (3 + 2n - P)y 2 - ty' 

6,1+1 - m=v ■ 

Substituting (20) and the expressions of b n ±i into (17), we indeed see that they 
satisfy the equation. Therefore, the system (6)-(7) is related to the Backhand 
transformation of Py. 

Finally we use Theorem 1.5 to obtain the parameters of Pm. Taking ki = 
1/2 from (21) we get the following values for different choices of the square 
roots: 

a = 2(1 + n - p), p = -2{n + p), (22) 
q = -2(3 + n - P), p = 2(7i + p), (23) 
a = 2(n - 1 + p), p = -2(2 + n-p), (24) 

a = -2(1 + n + P), P = 2(2 + n - p). (25) 

Taking fci = 1/2, a = t/2 we change variables t —> ^/2~i = z in (19) and take 
b n {a>) = B n (z). Then the functions 

- . . z + u(2P - 2n - 1 + zu) + zu' 
B n (z) = , 

fi / N 1 A n o a f 8(n+l) u' + u 2 + l 

5„(z) = - l + 2n-2# + z ; 

y ' A\ ' \z-u(zu + 2P-2n-b)-zu' u 

~ , . z + u(2p - 2n - 1 + zu) + zu' 
B n (z) = 

~ , , 1 /. „ no / 8(n+l) l + u 2 + u' 

BJz) = - l + 2n-2S + z ; r 

w 4V M V z + w(l + 2n + 2/?-zu)-zu' u 

give the transformations to the third Painleve equation for the function u = u(z) 
with parameters (22)-(25) respectively. 
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In case f3 = 1 we get from (22)-(25) either 

a = 2n, (j = -2(n+ 1) 

or 

a = -2(n + 2), /3 = 2(n+l). 

2.3 Initial conditions 

In this section we study the initial conditions (8), (10), (11) for the generalized 
Charlicr weight (4) on the lattice N, on the shifted lattice N + 1 — (3 and on 
the bi-latticc respectively. We show that the initial conditions correspond to 
classical solutions of Py- 

First we get that the equation (19) with n = has solutions satisfying 

a 2 b' Q + ab 2 Q - o(l - /3)b a - a 2 = 0, (26) 

where ' = d/da. The initial conditions (8), (10), (11) satisfy this equation. 
Moreover, by using the formulas from [1, 9], or computational software, it can 
be shown that the general solution of (26) is given by (11). We note that the 
Toda system (15) and discrete system system (6)— (7) are unchanged under the 
transformation (3 —> 2 — 0, a 2 = a 2 , b n — b n — 1 + (3. Also equation (26) is 
unchanged after this transformation with n = 0. 

Using the transformation (20) with ki — 1, a = t in (26) we get a first order 
second degree equation 

(y - 1)(1 + y(((3 - l) 2 y 2 - (At +((3- l) 2 )y - 1)) + 2t(y - l)y' - t 2 y' 2 = 0, 

for which the solutions satisfy Py with parameters A = (ft — l) 2 /2, B = 
— 1/2, C = 2, D = 0. Using the transformation y —> 1/y in Py we get 
equation (38.7) in [7] for special values of the parameters. Hence, the solutions 
of Py in terms of classical functions correspond to the initial conditions for the 
recurrence coefficients. 

Remark 4. The recurrence coefficients a 2 and b n can always be written [4] 
as ratio's of Hankcl determinants containing the moments of the orthogonality 
measure. However, the explicit determinant formulas of classical solutions for 
all the Painlcvc equations are known. 

2.4 Revisiting case (3 = 1 

It is shown in [3] that the monic polynomials orthogonal with respect to a 
semi-classical variation of the Charlier weight given by 

a k 

w(k) = a > 0, k e N, 
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have recurrence coefficients which satisfy Py with 

n 2 

A = -B = —, C = 0, D = -8. (27) 
8 

This weight can be obtained from weight (4) by taking j3 = 1. However, if wc 
take the parameter f3 = 1 in (21), wc get Py with 

(n + l) 2 

A = Q, B = - y 2 ' , C = 2h, D = Q, (28) 

which seems to be a contradiction. In the following we will show how to reduce 
(19) to Py with parameters (27). 

The following transformation exists between solutions of Py [7]. Let y = y(z) 
be a solution of Py with A = —B = a/A, C = 0, D = 4-f. Then the function 

m = ^r^, t = * 2 (29) 

Ay 

is a solution of Py with A = a, B = 0, C = 7, D = 0. So, if we start with 
a solution y of Py with parameters (28), we can first use the transformation 
y — > 1/y to get the parameters 

(n + l) 2 

A= y ' , B = 0, C = -2fci, Z? = 0, 
then the transformation (29) to get solutions with 

A=^t, B = -^t, C = 0, D = - 8kl . 
8 8 

Then taking fci = 1 and using Theorem 1.3 (e.g., a transformation from the 
example in case e = 1) we get a solution of Py with parameters (27). Explicitly 
this gives 

n + 7ny 2 - ny 3 - 2zy' - y(7n + 2zy') 2 

K[t) = m^) ■ 1 = z ' 

where y = y(z) is a solution of Py with (27). Equation (26) is also transformed 
to a first order second degree equation for which the solutions satisfy Py. Hence, 
we get that the solutions of Py in terms of classical special functions generate 
the initial conditions for the system (6)-(7). 



3 Discussion 



The recurrence coefficients for semi-classical weights are often related to the so- 
lutions of the Painleve type equations (e.g., [5, 6] and see also overview in [3]). 
The known results and our recent findings of the relations between the recur- 
rence coefficients of the (generalized) Charlier and Meixner polynomials and the 
fifth Painleve equation (12) with parameters (A, B,C, D) can be summarized 
as follows. We assume below that fcsN={0,l,...}. 
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1. (a) The weig ht a k /k\, a > 0, is a classical Charlicr weight and the recur- 

rence coefficients can be computed explicitly. 

(b) It is shown in [3] that the recurrence coefficients for the weight 
a k /(k\) 2 , a > 0, are related to the solutions of the fifth Painleve 
equation Py with parameters (n 2 /8, — n 2 /8, 0, — 8). Moreover, as 
shown in this paper, by using certain transformations, they are also 
related to solutions of Py with parameters (0, — (n + l) 2 /2, 2, 0) (this 
is the special case of 1(c)). 

(c) As shown in this paper, the recurrence coefficients for the weight 
a fc /((/3)fefc!), a > 0, (3 ^ 1, are related to solutions of Py with 
((/3-l) 2 /2,-(n+ l) 2 /2,2,0). 

2. (a) The weight ((3) k c k /k\, (3 > 0, < c < 1, is the classical Meixner 

weight and the recurrence coefficients are known explicitly. 

(b) It is shown in [3] that the recurrence coefficients for the weight 
(/3)fcC fe /(fc!) 2 , P > 0, c > 0, are related to solutions of Py with 
((/3-l) 2 /2,-(/3 + n) 2 /2,2n,-2). 

(c) It is shown in [6] that the recurrence coefficients for the weight 
(l)kC k /(fc!(/3)fc), c, /3,7 > 0, are related to solutions of Py with 
((7 - l) 2 /2, -(7 - /3 + n) 2 /2, ki(P + n), -fc 2 /2)), k x ^ 0. The case 
7 = 1 corresponds to classical Charlier weight on the shifted lattice 
N + 1 — (3. The case (3 = 7 gives the Charlier weight and the case 
f3 = 1 is given in 2(b). 

Therefore, we have essentially two different cases for Py: D = for the 
recurrence coefficients of the generalizations of the Charlier weight and D ^ 
for the recurrence coefficients of the generalizations of the Meixner weight. Note 
that the case D = can be reduced to the third Painleve equation [7]. 

It is interesting to study weights on shifted lattices. For instance, there is 
no shifted lattice for the weight a k /(k\) 2 since both lattices coincide. Since the 
lattices are determined by the zeros of the weight, one could have more than 
one lattice if one takes a fe /((/3)fc(7)fc(<5)fcfc!), which is zero for k = —1, k = —f3, 
k = —7 and k = —8. 

It is an interesting open problem to study recurrence coefficients for more 
complicated weig hts like a k /{k\f or {^) k (j3) k c k / {k\f and others. The Toda 
system will be the same but the discrete system will be much more complicated 
(maybe higher order discrete Painleve equations will be needed) . Will this lead 
to the Painleve equations (with more general parameters A, B, C, D depending 
on the parameters in the weight) or some higher order equations in some hier- 
archy of the Painleve equations? The direct calculations will be cumbersome 
and other techniques like using the Lax pairs or others will be needed to tackle 
such problems. 
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